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Abstract 

Various methods have been proposed to condition an electron 
beam in order to reduce its emittance effect and to improve the 
short-wavelength free electron laser (FEL) performance. In this 
paper, we show that beam conditioning does not result in a com- 
plete elimination of the emittance effect in an alternating-gradient 
focusing FEL undulator. Using a one-dimensional model and a 
three-dimensional simulation code, we derive a criteria for the 
emittance limitation of a perfectly conditioned beam that de- 
pends on the focusing structure. 
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1 Introduction 



A primary factor limiting the performance of short-wavelength, high-gain free 
electron lasers (FELs) is the electron beam transverse emittance. In order 
to generate transversely coherent radiation with the diffraction limited emit- 
tance Ar/47r (Xr being the radiation wavelength), the ideal electron transverse 
emittance e should be less than Xr/An for the efficient beam-radiation inter- 
action. This requirement is usually not met in the x-ray wavelength region. 
The spread in the transverse betatron motion degrades the resonant FEL 
interaction and limits the short- wavelength reaches of an x-ray FEL. 

To mitigate this problem, it has been proposed to "condition" an electron 
beam prior to the undulator by increasing each particle's energy in proportion 
to the square of its betatron amplitude [1]. This conditioning enhances the 
FEL gain by reducing the axial velocity spread within the electron beam 
generated over the undulator, due to both energy spread and finite transverse 
emittances. 

The original proposal [1] to condition the beam utilizes a set of cavities in 
TM210 mode immersed in a focusing lattice. Later approaches to the problem 
invoke a slow TM waveguide mode internal to the undulator [2] , and an en- 
ergy chirp in the beam in combination with a chromatic delay line [3, 4, 5, 6]. 
It has also been shown that an attempt to condition the beam on a short 
distance is accompanied by a head-tail focusing variation which can result in 
the large effective transverse emittance growth [4]. Although this emittance 
growth can be avoided [6], the pace of conditioning becomes much slower, 
and the required length of the conditioner considerably increases. Other 
proposals are based on using lasers and include Thomson backscattcring [7], 
and interaction of the laser with the electron beam in two dedicated undula- 
tors [8] . Recently, beam conditioning using nonlinear RF acceleration is also 
discussed [9]. 

In the FEL gain analysis of the original proposal [1], the undulator nat- 
ural focusing is assumed to confine the electron beam in both transverse 
directions. This results in the complete elimination of the emittance effect 
for a perfectly conditioned beam. In this paper, we study a more practical 
situation when the transverse focusing is provided by alternating-gradient 
quadrupoles in the undulator, as found in typical x-ray FEL designs. We 
show that although conditioning eliminates the effect of the average slippage 
of electrons relative to the radiation phase, there remain phase oscillations 
with the period given by the focusing lattice. These phase oscillations can 



affect the FEL performance and impose a limitation on the beam emittance 
even for a perfectly conditioned beam, albeit this emittance limitation is 
much relaxed as compared to a unconditioned beam. 



2 Emittance effect in a strong focusing undu- 



Although a magnetic undulator can provide natural focusing in both trans- 
verse planes [10], the focusing strength is typically too weak for the high- 
energy electron beam that drives an x-ray FEL. Thus, alternating-gradient 
quadrupole magnets are inserted in undulator breaks to provide the necessary 
strong focusing, usually in the form of a FODO lattice (consisting of repeti- 
tive focusing-undulator-defocusing- undulator cells). The horizontal betatron 
motion is given by 



where is the horizontal action of the electron, and (3^ are the Twiss 
parameters, and denotes the angle of the orbit with the z axis. The second 
term in Eq. (1) is ignored in Ref. [11] under the smooth approximation. Its 
importance in FEL dynamics and emittance compensation (i.e., conditioning) 
is pointed out in Ref. [12]. 

To avoid a large beam size variation in the undulator, the FODO lattice 
is usually designed to have a small phase advance per cell. In this case, the 
average betatron wavelength 27r/9 is much larger than the FODO cell length 
Lc- In a first approximation the beta function is constant over the length of 
the undulator, and the accumulated betatron phase is 



where (f)x is the initial electron phase. More importantly for the emittance 
compensation, such a FODO lattice has [12] 
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Figure 1: Variation of the beta function along the distance of the lattice for 
the phase advance per cell of 15 degrees. The derivative d/S^/dz is close to 
the values ±2, but the deviation of (3^ from the average value /3 is relatively 
small. 

with the sign alternating per half cell. Fig. 1 illustrates the variation of the 
beta function in a FODO lattice with 15 degrees phase advance per cell. 

Under this short-cell-length approximation, the square of the horizontal 
divergence is 



Here the signs ± becomes =F due to the focusing-defocusing asymmetry. We 
have also assumed that the average beta function is the same in both direc- 
tions. 

The FEL interaction is a resonant phenomenon that depends critically on 
the evolution of the electron phase relative to the co-propagating radiation 
field. The so-called ponderomotive phase is defined as 6^ = {kr + ku)z — 
ckri, where kr = 27i/Xr is the radiation wavenumber, fc„ = 2n/Xu, A„ is the 
undulator period, and i is the electron arrival time at the location z averaged 
over the undulator period. The rate of the phase change is 





Similarly in the vertical y direction, we have 
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where is the undulator-period-averaged longitudinal velocity and is given 

by 

Here K is the undulator parameter and is approximately constant. Its weak 
dependence on x and y can be neglected as long as the external focusing of the 
FODO lattice is much stronger than the natural focusing of the undulator. 
Making use of the resonant condition 

we can write Eq. (6) as 
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Since the last term of Eq. (9) is oscillatory with the FODO lattice period 
(cell length), the main accumulating effects on the ponderomotive phase are 
the first two terms. Neglecting this oscillatory term, the phase equation in 
such a FODO cell is identical to that in a natural focusing undulator [13], as 
noted in Refs. [12, 14]. For a beam with finite energy spread and emittance, 
the first two terms introduce phase slippage of the electron relative to the 
radiation and result in the phase spread of the beam. Denoting () as the 
average over the beam, we have {J^ = Sx and {Jy) = Sy, where ex,y = £n/lr 
are the beam transverse emittances. In order to not significantly degrade the 
EEL performance, the emittance-induced phase spread over one EEL power 
gain length Lq should be less than unity, i.e.. 

This situation changes if the initial electron energy can be conditioned to 
its transverse betatron amplitudes such that 

2k„^ = k^MM, (11) 

It P 



the dominant emittance effect (i.e., the second term in Eq. (9)) is then re- 
moved from Eq. (9). Averaging over and Jy^ the (transversely) correlated 
energy spread required to condition a beam with the normalized emittance 

(A.o).^^|. (12) 

However, when the beam is perfectly conditioned to satisfied Eq. (12), 
the last oscillatory term in Eq. (9) is no longer negligible and can play a 
limiting role with a large enough emittance. We discuss quantitatively this 
remaining emittance effect on the FEL performance for a conditioned beam 
in the following two sections. 

3 ID model of the ponderomotive phase os- 
cillation 

The ponderomotive phase oscillations in Eq. (9) depend on the transverse 
variables through J^^ Jy, d>x and 0,^ and is a three-dimensional (3D) problem. 
To isolate the oscillation effect and to simplify the problem, we study a 
heuristic one-dimensional (ID) model with the following phase equation 

f^2kj+^fiz), (13) 

where S = {Aj — Ajq) / jr is the FEL- induced energy change, and f{z) repre- 
sents the oscillatory behavior introduced by a conditioned beam in the FODO 
lattice (with a cell length much smaller than the betatron wavelength), i.e.. 



f(^)=)+'^ when (n- l)Le < ^ < (2n- 1)Y, 

' -1 when (2n- 1)^ < ^ < nLe, ^ ^ 



and n = 1,2,3.... In the absence of the FEL interaction (i.e., when 6 = 0), 
the phase oscillates between and with the maximum phase deviation 

Using the FEL Pierce parameter p [15], we introduce a scaled distance 
T = 2kupz and a scaled energy -q — 5/ p. Equation (13) becomes 

f = r; + e, with e = ^^/(r). (16) 



The FEL-induced energy change is 

— — ae*^ + complex conjugate , (17) 
or 

where a is the slowly-varying radiation field amplitude (properly scaled by 
p) . Neglecting any transverse dependence and considering the radiation field 
at the resonant frequency ckr, the ID field equation is 

| = -(e-^). (18) 

This set of coupled equation can be solved by averaging over the fast 
oscillation if its period is much smaller than the field gain length, similar 
to the undulator-period averaging procedure for a planar undulator. To 
illustrate this process, we define 



^ = ^-y ar')dr', (19) 
Equation (16) can be written as 

Equations (17) and (18) become 

^ =ae'^exp (i [ ^(T')dT') + complex conjugate, (21) 
dT \ Jo J 

^ = -(e-Vxp(-^ [ i{T')dT^ . (22) 



\ -'0 / 

Treating r/, ^, and a as slowly-varying variables, we can average Eqs. (21) 
and (22) over the oscillation period to obtain 

— —ae^^A{6Q) + complex conjugate, (23) 
or 

^ = - (e-^V(-^o) , (24) 
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Figure 2: ID FEL power growth rate under a rapid phase oscillation with 
the maximum phase deviation 6q. Solid curve is the solution of Eq. 26, and 
the symbols represent ID simulation results. 
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Equations (20), (23) and (24) are the FODO-cell averaged FEL equations 
and can be solved with the usual techniques (see, e.g., Ref. [15]). Assuming 
that a(r) oc e"*^"^, we obtain a cubic equation for the complex growth rate 
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The FEL power growth rate 2Im/x = v^|^(^o) P''^ versus the maximum phase 
deviation is shown in Fig. 2. It agrees well with the ID simulation results 
by solving the original set of Eqs. (16), (17) and (18) for an electron beam 
without any initial energy spread. The FEL growth rate is degraded when 



h ^ T 
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Table 1: GENESIS PEL simulation parameters of conditioned beams in a 
strong focusing uudulator. 



Parameter 


Symbol 


Value 


electron energy 




20 GeV 


flattop bunch current 


Ipk 


3 kA 


transverse norm, emittance 




0.1 to 19 /xm 


conditioned energy spread 




given by Eq. (12) 


uncorrelated energy spread 







undulator period 




3 cm 


undulator parameter 


K 


4.67 


average beta function 


P 


4.8/9.4/19 m 


FODO cell length 


Lc 


2.64/5.04/9.83 m 


PEL wavelength 


Ay. 


1.0 A 



4 3D GENESIS simulations 

The above ID model shows that when the maximum phase oscillation ampli- 
tude exceeds unity, the PEL gain will be degraded as compared to the ideal 

case. Therefore, for a perfectly conditioned beam that satisfies Eq. (12) in a 
PODO lattice, the PEL performance will still be affected by the emittance- 
induced phase oscillations when the normalized emittance exceeds a crit- 
ical emittance given by 



Since this emittance criteria is derived with a heuristic ID model, we 
examine it using the 3D PEL code GENESIS [16]. The electron beam and 
undulator parameters are given in Table 1. A typical PEL power evolution 
(P versus z) is shown in Pig. 3 (blue dashed curve). The local power growth 
rate l/P{dP/dz) is also plotted in the same figure (green sohd curve) and is 
oscillatory with a period equaling to Lc/2 as a result of emittance- induced 
phase oscillations. 

We extract the PEL power gain length by averaging a relatively constant 
l/P{dP/dz) over many oscillation periods in the GENESIS simulations. The 
resulting gain length versus the normalized emittance is plotted in Pig. 4 for 
;9 = 19 m and for three different PODO cell lengths. We also compare the 
simulation results with the theoretical gain length for a conditioned beam 
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Figure 3: (Color) GENESIS amplifier run of a conditioned beam for L^. = 
5.04 m, /? = 19 m and En = 10 /im. The blue dashed curve is the PEL 
power P at Ar =1 A, while the green solid curve is the local growth rate 
l/P{dP/dz). 

in the absence of any phase oscillation. The latter is obtained by solving 
the usual 3D eigenmode equation (see, e.g., Ref [17]) for the fundamental 
Gaussian mode without any energy and angular spreads. As shown in Fig. 4, 
the simulation results agree with the 3D theory that does not take into ac- 
count phase oscillations up to the critical emittances predicted by Eq. (28). 
For emittances exceed these critical values, the gain lengths extracted from 
simulations start to increase faster than the ideal case. 

We have also performed GENESIS simulations by varying the average 
beta function (3 while keeping the same focusing structure (with Lc = 2.64 m). 
Figures 5 and 6 show the gain length comparisons for Lc = 2.64 m at /? = 
9.4 m and 4.8 m, respectively. Combining with Fig. 4, we see that Eq. (28) 
is reasonably accurate in predicting the critical emittances for different beta 
functions and FODO lattices. 

As the ratio of the average beta function to the cell length becomes too 
large (e.g., when /3 = 39 m and = 2.64 m or when /5 = 19 m and 
Lc = 1.44 m), the gain length extracted from the GENESIS simulation starts 
to deviate from the 3D theory (without any phase oscillation) at a smaller 
emittance than that predicted by Eq. (28). A plausible explanation is given 
as follows. In x-ray FELs the transverse coherence builds up slowly as many 
higher-order transverse modes are excited by an electron beam with a rela- 
tively large emittance. To some degrees a periodic modulation in the electron 
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Figure 4: (Color) FEL power gain length from GENESIS simulations at 
/5 = 19 m for Lf. = 2.64 m (blue solid curve 1), = 5.04 m (green solid 
curve 2), = 9.83 m (red solid curve 3), and from 3D theory without any 
phase oscillation (black dashed curve). The colored symbols represent the 
corresponding critical emittance for each lattice (determined by Eq. (28)) 
when the phase oscillation is expected to increase the gain length. 
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Figure 5: (Color) FEL power gain length from GENESIS simulations at 
f3 = 9.4 m for Lc = 2.64 m (red solid curve), and from 3D theory with- 
out any phase oscillation (black dashed curve). The symbol represents the 
corresponding critical emittance (determined by Eq. (28)) when the phase 
oscillation is expected to increase the gain length. 
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Figure 6: (Color) FEL power gain length from GENESIS simulations at 
jS = 4.8 m for = 2.64 m (red solid curve), and from 3D theory with- 
out any phase oscillation (black dashed curve). The symbol represents the 
corresponding critical emittance (determined by Eq. (28)) when the phase 
oscillation is expected to increase the gain length. 

beam size along the undulator is beneficial to the build-up of the transverse 
coherence (namely higher order modes are suppressed by a varying beam en- 
velope [18]). When the ratio of the beta function to the cell length becomes 
too large, this higher-order suppression mechanism is not effective since the 
beam size is almost constant, and the growth of the radiation is contributed 
from both the fundamental mode as well as higher-order modes. In this case, 
the comparison between simulations (including higher-order modes) and the 
theory (for the fundamental mode) is no longer valid. 

5 Summary 

In this paper, we study the FEL performance for a perfectly conditioned 
beam in a strong focusing undulator consisting of FODO cells. We develop 
a heuristic ID model for the emittance-induced phase oscillation and show 
that the maximum phase deviation should be less than unity in order not 
to degrade the FEL gain. This criteria limits the maximum emittance that 
may be conditioned in such a focusing channel and is confirmed by the 3D 
GENESIS simulations over a reasonably wide parameter range. Therefore, 



instead of the usual emittance criteria for the unconditioned beam as given 
by Eq. (10), the emittance of the conditioned beam in a FODO lattice with 
an average beta function j3 and a cell length Lc must satisfy 



That is, the emittance requirement for a conditioned beam is much relaxed 
but not without any limitation. 



6 Acknowledgments 

We thank A. Zholents for useful discussions. This work was supported by 
Department of Energy contract DE-AC02-76SF00515. 



References 

[1] A. Sessler et al, Phys. Rev. Lett., 68, 309 (1992). 

[2] P. Sprangle et a/., Nucl. lustrum. Methods A 331 (1993) 6. 

[3] N. A. Vinokurov, Nucl. Instrum. Methods A 375 (1996) 264. 

[4] P. Emma and G. Stupakov, Phys. Rev. ST Accel. Beams, 6, 030701 
(2003). 

[5] P. Emma and G. Stupakov G. Proceedings of 9th European Particle 
Accelerator Conference, p. 503 (2004) 

[6] A. Wolski, G. Penn, A. Sessler, and J. Wurtele Phys. Rev. ST Accel. 
Beams 7, 080701 (2004) 

[7] C. B. Schroeder, E. Esarey, and W. P. Leemans Phys. Rev. Lett. 93, 
194801 (2004) 

[8] A. Zholents, Phys. Rev. ST Accel. Beams 8, 050701 (2005) 

[9] G. Stupakov and Z. Huang, in Proc. of the PEL '05, Stanford, CA, USA, 
p. 176 (2005). 



[10] E. Scharlemann, J. Appl. Phys. 58, 2154 (1985). 

[11] L.-H. Yu, C. Huang, D. Li, and S. Krinsky, Phys. Rev. A 51, 813 (1995). 

[12] S. Reiche, Nucl. Instrum. Methods A 445, 90 (2000). 

[13] In a natural focusing undulator, ax,y = in Eq. (1), but the trans- 
verse variations of the undulator parameter K together with the first 
term in the bracket of Eq. (1) result in the phase equation d9/dz = 
2kuAj/'jr — kr{Jx + Jy)/ l^m whcrc /9„ = 7rA„/(7ri^r) is the natural focus- 
ing beta function. See also Ref. [10]. 

[14] E. Saldin, E. Schneidmiller, and M. Yurkov, Nucl. Instrum. Methods A 
475, 86 (2001). 

[15] R. Bonifacio, C. Pellegrini, and L. Narducci, Opt. Commun. 50, 373 
(1984). 

[16] S. Reiche, Nucl. Instrum. Methods A 429, 243 (1999). 
[17] M. Xie, Nucl. Instrum. Methods A 445, 59 (2000). 
[18] S. Reiche, Ph. D. dissertation, DESY (1999). 



